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Abstract 

It is of interest in a variety of contexts, and in particular in the arrival time problem, to consider 
the quantum state obtained through unitary evolution of an initial state regularly interspersed with 
periodic projections onto the positive x-axis (pulsed measurements). Echanobe, del Campo and 
Muga have given a compelling but heuristic argument that the state thus obtained is approximately 
equivalent to the state obtained by evolving in the presence of a certain complex potential of step- 
function form. In this paper, with the help of the path decomposition expansion of the associated 
propagators, we give a detailed derivation of this approximate equivalence. The propagator for the 
complex potential is known so the bulk of the derivation consists of an approximate evaluation of the 
propagator for the free particle interspersed with periodic position projections. This approximate 
equivalence may be used to show that to produce significant reflection, the projections must act at 
time spacing less than h/ E, where E is the energy scale of the initial state. 

PACS numbers: 
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I. INTRODUCTION 



A number of physically interesting situations in quantum theory concern the question of 
what happens to an initial quantum state l'^) which is acted on by a sequence of projection 
operators interspersed with unitary evolution for time e, 



(Here, there are n projection operators and t — {n + l)e and we use units in which h—1). 
Such an object describes pulsed measurements but it also crops up in the decoherent histories 
approach to quantum theory, where the right-hand side is the amplitude for a quantum- 
mechanical history [1-6]. 

We are interested in the specific case of a free particle with P taken to be the projection 
onto the positive x-axis, P — 6{x). For sufficiently small e, Eq.(l.l) is then a candidate for 
the amplitude to remain in x > during the time interval [0, r], an object that is of interest 
in the arrival time problem [7-10]. 

It is of interest to explore the properties of this amplitude for a range of values of the 
time spacing e. It is known that as e ^ 0, wc approach the Zeno limit, in which the state 
becomes entirely confined to the Hilbert subspace of states with support only in a; > 0, so 
that an incoming wave packet from the right is totally reflected [11-13]. However, it is of 
greater physical interest to explore the regime of non-zero e, in which the system is monitored 
sufficiently well to get some idea of whether the particle is in x > 0, yet not monitored so 
much that an incoming state is significantly refiected at x — 0. An important question in 
this regime is to determine the value of e for which refiection becomes significant. For an 
initial state with energy width AH, a timescale held to be significant is the Zeno time. 



which is the timescale on which the state becomes significantly different from its initial 
value under unitary evolution [11]. For a wave packet of momentum p and spatial width a, 
the Zeno time is of order ma/p which is the timescale on which the wave packet crosses the 
origin. This indicates that the Zeno time for wave packets is an essentially classical timescale 
and, in Eq.(l.l), relates only to the rate of removal of probability through projection. By 
contrast, refiection in Eq.(l.l) arises as a result of the increase in uncertainty in momentum 



\Mr)) 



= e-'^'Pe-'^'-.-Pe-'^'lip) 



(1.1) 
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resulting from position projection, an obviously quantum process, so one would expect it to 
have a different timescale, which could be much shorter than then Zeno time. It would be 
of interest to compute this timescale. One reason it is important is that there appears to be 
interesting physics very close to the Zeno limit [14-16]. 

A significant result in this area is due to Echanobc, del Campo and Muga, who claimed 
that for finite e the string of projection operators in Eq.(l.l) is approximately equivalent to 
evolution in the presence of a complex potential [14]. That is, 

^-tHep^-tm . . . p^-tm ^ _ Voe{-x)r) (1.3) 

They assert that this result is valid if, for a given e, Vq is chosen such that two conditions 

Voe > 1 (1.4) 
Vo > AH (1.5) 

are satisfied. This is a very useful result since Eq.(l.l) is not easy to evaluate analyti- 
cally but the Schrodinger equation with a complex step potential in Eq.(1.3) can be solved 
straightforwardly. Furthermore, such complex potentials have been studied extensively in 
the literature and can often be linked to particular detection methods [17]. 

Given the connection Eq.(1.3), one can determine the conditions under which reflection 
becomes important. Known results on scattering with the complex potential in Eq.(1.3) 
show that, for an incoming state with energy scale E, reflection becomes signiflcant when 
Vq > E [8, 9]. Reflection is avoided, therefore, when Vq -C E, or equivalently, from Eq.(1.4), 
when 

e » 1 (1.6) 

This is much less than the Zeno time for a state strongly peaked in energy. There is therefore 
an interesting regime, namely 

1 « ^ « ^ 

in which the projections in Eq.(l.l) are sufficiently frequent to have a significant effect on 
the system, yet not so frequent that there is significant reflection. 

Eq.(1.6) is a very useful result, but it has been derived on the basis of the claimed 
approximate relationship Eq.(1.3). The derivation of Eq.(1.3) given by Echanobe et al 
is very plausible (and was also hinted at by AUcock [8]), but it is rather heuristic, and the 
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deduced connection Eq.(1.4) between Vo and e is rather loose. There is therefore considerable 
scope for a more detailed and substantial derivation. 

The purpose of this paper is to give a more substantial derivation of the equivalence 
Eq.(1.3) and to deduce a more precise relationship between Vq and e. We will do this by 
computing the configuration space propagators associated with each side of Eq.(1.3) and 
show that they are approximately equal in certain regimes. The propagator associated 
with the complex potential is in fact known already, so the bulk of the work consists of an 
approximate evaluation of the propagator associated with pulsed measurements, the left- 
hand side of Eq.(1.3). 

This is certainly not a rigorous mathematical proof of Eq.(1.3), involving operator norms, 
error bounds and the like, but a theoretical physicists style of proof involving the approximate 
evaluation of propagators. A rigorous proof would certainly be of interest to construct and 
the work described in this paper may give some hints in that direction. 

We begin in Section 2 with a brief summary of the derivation of Echanobe et al, with a 
small extension of it that turns out to be important and yields an equality relating Vq and 
e, thereby improving on Eq.(1.4). We then in Section 3 give a detailed formulation of the 
problem and how we solve it. The key idea is to use the path decomposition expansion in 
which the propagators associated with each side of Eq.(1.3) are factored across the surface 
X — [18-20]. The problem of proving the equivalence Eq.(1.3) thereby reduces to proving 
it for propagation between points ai x — for different times. Since the propagator for 
the complex potential is known, the main work is to evaluate this propagator for pulsed 
measurements. This is actually rather difficult to do directly, but a good approximate 
analytic expression can be obtained for it by attacking the problem from a number of different 
angles. This is described in Sections 4, 5, 6 and 7. In Section 8 we give a detailed discussion 
of the timescales involved for the approximations to be valid. We summarize and conclude 
in Section 9. 

II. REVIEW AND EXTENSION OF EARLIER WORK 

We first review and extend the derivation of Echanobe et al [14]. They first note that 




(2.1) 
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where, recall P — 9{x) and P — 1 — P. It follows that 

P = e{x) ^ exp {-Voe{-x)e) (2.2) 

as long as the parameter 



a = Voe (2.3) 



is sufficiently large that 



e-" < 1 (2.4) 
Eq.(1.3) then follows from the approximate equivalence, 

exp {-iHe) exp {-Voe{~x)e) ^ exp {-iHe - Vo9{-x)e) (2.5) 

which will hold as long as 

Voe'\{[H,e{-x)])\^i (2.6) 

where the average is taken in the initial state. Echanobe et al put an upper bound on the 
left-hand side using the Schrodinger- Robertson inequality and Eq.(2.6) may then be written 
in either of the two equivalent forms 

a^AH < Vo (2.7) 

or 

ae « ^ (2.8) 

which implies that the time between projections is much less than the Zeno time, the typical 
timescale on which the state undergoes significant change [11]. (The conditions Eqs.(2.3), 
(2.4), (2.7) are a more precise version of the originally stated conditions Eqs.(1.4), (1.5)). 

Since the parameter a need only satisfy an inequality, Eq.(2.4), the relationship Eq.(2.3) 
between e and Vq is not uniquely determined. However, it turns out that the above derivation 
can be extended somewhat to give an equality between e and Vq. This turns out to be relevant 
to the more substantial derivation given in the rest of this paper. 

The above result may be written, 

exp {-iHe - Vo9{-x)e) ^ exp {-iHe) P (2.9) 
under the conditions given above. Now suppose a is an integer and write e = ae', where 

Voe' = 1 (2.10) 



Then, since P — P'^,we may approximate Eq.(2.9) by 

eicp{-iHe-Voe{-x)e) ^ e-'^'' Pe-'"'' P ■ ■ ■ e''"'' P 

= (^e-'"''py (2.11) 

as long as the contribution from the commutator terms between e~'^^' and P is sufficiently 
small. There will be of order such terms, hence the error in this approximation is of order 
a^e'AH and from Eq.(2.7), this error is much less than 1. The key point here is that even 
the longer timcscalc e is still much less than the Zeno time and, since nothing changes on 
this timescale, there is essentially no difference between Eqs.(2.9) and (2.11). 

We therefore see that the desired result Eq.(1.3) actually holds for much smaller time 
steps, defined by the equality Eq.(2.10). This is different to the original claim of Echanobe 
et al, since they require the inequality Eq.(1.4). However, the above argument shows that 
the restriction Eq.(1.4) is in fact stronger than necessary and in the following pages our 
more detailed derivation will show that Eq.(1.3) does indeed hold with a timespacing of 
order I/Vq. 

III. DETAILED FORMULATION OF THE PROBLEM 

In this paper, we will prove the relationship Eq.(1.3) in a much more substantial way by 
proving the approximate equivalence of the propagators 

gv{xi,r\xo,0) = {xi\ exjp {-iHr - Vo9{-x)r) \xo) (3.1) 
gp{xi,T\xo,0) = {xi\e-'"'"Pe-'"'...e-'"'Pe-'"'°\xo) (3.2) 

for some relationship between the parameters e and Vq, to be determined. Note that in 
Eq.(3.2), we have chosen the initial and final time spacings to be eo and e„, with 

T = (n - l)e + eo + e„ (3.3) 

This turns out to be necessary for the most general proof of Eq.(1.3). For the special case 
en = e = eo, we may also write 

gp{xi,T\xo,0) = (xi|e-^^^P(ne)---P(2e)P(e)|xo) (3.4) 
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Each of the above propagators may be represented by a path integral, 

gv{xi,T\xQ,0) — J Vxexp^ J dt -mx'^ + iVo9{—x) ^ (3.5) 

gp{xi,T\xo,Q) = J VxexTp^ J dt^mx^^ (3.6) 

where in both cases the paths are from x{0) — xq to x{t) — xi and in the second case, 
Eq.(3.6), P denotes that the paths are restricted to be in the positive x-axis at times 
t = eo + {k — l)e, k = 1,2 ■ ■ ■ n. 

A closely related object that will be important is the restricted propagator, 

gr{xi,T\xQ,0) = j Vxex])(i f dt-mx^j (3.7) 
Jx{t)>o V 2 / 

where again the paths are from ,x(0) = Xq to x(r) = Xi but with x(t) > for all times 
in [0,r]. This is clearly equivalent to gp in Eq.(3.6), in the limit n — )■ oo, e — ?> with r 
constant. If we take the same limit in the equivalent expression for gp, Eq.(3.2), one obtains 
the following convenient operator form of the restricted propagator: 

gr{T, 0) = P exp{-iPHPT) (3.8) 

The restricted propagator satisfies the Schrodingcr equation in a; > subject to the boundary 
conditions that it vanishes when either end of the propagator sits on a; = 0. For the free 
particle, considered here, one can easily solve for the restricted propagator using the method 
of images and the result is 



/ ffi \ 1/2 

gr{xi,T\xo,0) = ['^~) 0{xi)e{xo) 



X 



im(xi — xoY\ fimixi+xoY 



exp — exp 



2t / "V 2t 



(3.9) 



The restricted propagator, in any of the above forms, describes the regime of "Zeno dynam- 
ics" , in which all states are confined entirely to the Hilbert subspace of states with support 
only in X > [13]. 

The restricted propagator plays an important role here since not only does gp ^ gr 'nil 
the limit n — )■ oo, e — )■ with r constant, but also g^y — )■ as Vq — >■ oo. It follows that 
gv and gp become arbitrarily close to each other for sufficiently large Vq and n, since they 
both tend to the same limit. This is the underlying reason why we expect the approximate 
equivalence Eq.(1.3) should hold. 



The propagators gp and gv may be decomposed using the path decomposition expansion 
(PDX), in which a propagator, regarded as a sum over paths, is spht into propagation 
corresponding to paths (or sections of paths) entirely in x > 0, paths entirely in x < and 
paths starting and ending on the boundary [18-20]. There are two cases to consider. 

We consider first the case in which the initial and final points are on the same side of 
the surface in x > 0. The set of paths from initial to final point may be partitioned into 
paths that never cross the origin, represented by the restricted propagator gi^, and paths 
that always cross. The paths that cross have a first crossing at ti and a last crossing at ^2- 
The corresponding PDX for any propagator g has the form 

1 r Z"*^ dg dg 

+ gr{xx,T\xo,Q) (3.10) 

This is depicted in Figure 1. The first term in this expression, describing paths which cross 
the origin, contains the derivative of gr because this corresponds to paths which never cross 
X — Q but end or begin at x = (recalling that gr itself vanishes if either end is on x = 0) 
[18]. 

Eq.(3.10) holds for both gp and gy^ but the restricted part g.^. is the same in each case 
and equal to Eq.(3.9) above, since the restriction on paths defining gp and the presence of 
the complex potential in gy in Eqs.(3.5), (3.6) are both redundant if x{t) > 0. Eq.(3.10) 
may therefore be simplified using 

^(a:,ti|a:o,0)U = 2^(0,ti|a:o,0) (3.11) 

which follows from Eq.(3.9), where gf denotes the free particle propagator. 

It follows from the above that gp and gv could differ only in terms of their propagation 
along a; = 0, hence to prove the approximate equivalence of the propagators Eqs.(3.1), (3.2), 
we need to prove that 

gp{^M^M) ^ gvi^M^M) (3.12) 

Note that, unlike 5iv(0, ^2(0, ti), 5'p(0, t2|0, ti) is not covariant under time translation, that 
is, 

^p(0,t2|0,ti) 7^^p(0,t2-ti|0,0) (3.13) 

although it has an approximate covariance on timescales much greater than e, as we will see 
below. 



8 



The second case is that in which the initial and final points are on opposite sides of the 
surface, xi < and xq > 0. In this case we have for gv, 

1 Z"*^ da 

gv{x,,r\xo,0) = ^ dh e'^^^-"*^) dt, -^{x,,r\x,h)l^^ 

X gv{0,t2\0,h)^{x,h\xo,0)l^^ (3.14) 

since the last section of the paths is in a; < where the complex potential acts. For gp, if 

Xi < 0, since the paths must be in x > at the given discrete set of times, the last crossing 
time t2 cannot be less than the last time r — e„ at which the projectors act, hence 

rt2 



1 r n dg-r 

gp{xi,T\xo,0) = J dt2 J dti -^{xi,T\x,t2 



x=Q 



X gp{0,t2\0,ti) -^{x,h\xo,0)l^^ (3.15) 

Again we will have to prove that Eq.(3.12) holds, but these two expressions Eqs.(3.14), 
(3.15) also differ in the form of the t2 integral. However, since Vo€ fa 1, the form of the 
exponential in Eq.(3.14) effectively squeezes t2 to lie approximately within e of r, so we have 



/ dt2 e-^o^^-*^) « / 

Jo J T—en 



dt2 (3.16) 



and Eqs.(3.14) and (3.15) are approximately the same. 

We see that both cases reduce to proving Eq.(3.12). The propagator along the boundary 
for the complex potential is known [21], and is given by 

The main purpose of the remainder of this paper is to calculate the propagator with projec- 
tion operators along the boundary, 

gp{0,r\0,0) = {0\e-'"'-Pe-'"' ■ ■ ■ e-'"'Pe-'^'°\0) (3.18) 

and show that the approximation Eq.(3.12) holds, under conditions to be determined. (Here, 
|0) denotes a position eigenstate \x) at x = 0). 

In what follows, our main result is to show that, to a good approximation, 

^p(0,i|0,0)«. (^) fp{t) (3.19) 
9 



where fp{t) is a kind of saw-tooth function - a piecewise hnear function with peaks imme- 
diately followed by troughs at tk — eo + {k — l)e, k — 1,2, ■ • ■ . Approaching tk from below, 
there is a peak of value 

and approaching tk from above there is a trough of half that size. That is 

/pW = ^|^^ + %^ for tk-^<t<tk, k = 2,3... (3.21) 

and 

fp{t) = 1 for < t < eo (3.22) 

The functions fp{t) and fv{t) are shown in Figure 2. We see that fp{t) oscillates with period 
e about fv{t), as long as we choose Vo so that fv{t) lies between the peaks and troughs of 
fp{t). That is, for large k, we require that 

111 

— < — < - (3.23) 

Since t ~ ke, fvif) will lie approximately midway between the peaks and troughs of fp{t) if 

Voe ^ \ (3.24) 

Recalling that the propagator is attached through the PDX Eq.(3.10) to an initial state, 
the oscillations, and hence the differences between gp and gv, will be smoothed out as long 
as e is chosen to be smaller than the timescale of variation of the initial state. With some 
qualifications (discussed further in Section 8), this timescale is the Zeno time, 1/AH. The 
desired approximation Eq.(3.12) will therefore hold in a time-averaged sense, and we thus 
have significant agreement with the extended version of the original argument of Echanobe 
et al described in Section 2. 

In the following sections, we evaluate Eq.(3.18) and confirm the form Eq.(3.21) of fp{t). 
In Section 4 we evaluate Eq.(3.18) exactly for the cases of one, two and three projections. 
We also show why in general the troughs of fp{t) are exactly half the size of the peaks 
by considering the limit e„ — ?■ of Eq.(3.18). In Section 5, we use a lattice method to 
derive the magnitude of the peaks of fp{t) for large k. These results are substantiated in 
Section 6, where we use an S-matrix expansion to derive some exact results for a certain 
time-averaged version of Eq.(3.18). In Section 7, we fill in some of the gaps in these regimes 
and approximations by computing Eq.(3.18) using numerical methods. 
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We will find in the numerical and analytic solutions that the function interpolating be- 
tween the troughs and peaks of fp{t) is not in fact a linear function in general. It is inter- 
esting however, that, because of the slow time variation of the initial state in comparison 
to e, the precise form of this function turns out to be unimportant, and this is what makes 
the problem more tractable than one might expect. All that is important is the location of 
the peaks and troughs of the saw-tooth function and the period of their oscillation. It is 
this separation of timescales that also restores an approximate time translation covariance 
to gp, even though it does not hold exactly in general, Eq.(3.13). 



IV. EXACT ANALYTIC RESULTS 



In this section we carry out an exact evaluation of Eq.(3.18) for the case of one, two and 
three projections. For simplicity, we take the initial time interval eo to be e and for the case 
of three projections we are able to carry out the calculation only in when the final time 
e„ = e. We thus compute the objects 



gp{Q,t\Q,Q)^ { 



(0|e-'^*|0), 

(0 1 e-^^(*-2^) p^-iHep^-iHe | ^ 



if < i < e 
if e < t < 2e 
if 2e < t < 3e 



(4.1) 



{0\e-iHep^-iHep^-iHep^-iHe^ ^ if t = 4e 



We will show that 



^p(0,t|0,0) 



(2^)' ifO<t<6 
H^)' iie<t<2e 
i (1 + f ArcTan[(i - 2e)/i]^) (^)^ if 2e < i < 3e 



(4.2) 



1 { m 
4 



thereby confirming the approximate form of the saw-tooth function, Eq.(3.21), for the first 
few projections. These expressions all have the property that the value of gp{0, t\0, 0) drops 
to half its value immediately after a projection and we show why this is true in general at 
the end of this section. 

To compute Eq.(4.1) it turns out to be useful to consider more general objects, in which 
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the projections may be P, P or the identity. We therefore consider the object 

Here C = {+,—,0...} symbohcally stands for the integration ranges of the Xk, eg C — 
{+, — ,0...} means < xi < oo, — oo < X2 < 0, — oo < X3 < oo etc. Performing a Wick 
rotation, e — > —ie and changing variables yields 

dyi...dyn 



/ jji \ 1/2 c 
<,P,c(0,-«|0.0) = (-) 



c (ire)»/2(t-Ti£)i/2 

:= (4.6) 

The cases of no projection and one projection are trivially evaluated, and we easily obtain 
the first two equations in Eq.(4.2). For the case of two projections, we need to evaluate the 
object 

r++ ei,e2,e3 = / — ^=exp (4.6) 

Jo 71^616263 V ei 62 €3/ 

(a slight generalization of the Tc defined in Eq.(4.5)). Changing variables to r — i/i, y — 
Z/2/yi gives, 

r++ ei,e2,e3 = / exp -r^ <^ - + ^ ^ + ^ 

dy_ 1 

27r (€263 + 6163(1 - y)2 + 61622/^) 
/o 27r a + % + 

Where a = 63(61 + 62), h — —26163, c = 61(62 + 63). Now use the change of variables. 



(4.7) 



u^^^ + ^cy (4.8) 



to obtain 

1 



r++(6i,62,63) = — ./-f- / du (4.9) 
27r V £1630 ^2 + 4ac_0_ 

' 2vc 4c 
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Noting that 4ac — 6^ > 0, and using the standard integral [22] 



/ 



du ^ = -ArcTan (-) (4.10) 



we find finally 



r++(ei,e2,e3) = ^ = tt + 2 ArcTan ( ./ (4.11) 

47rv/(ei + £2 + £3) V VV ^2(61 + £2 + £3)/ / 

Setting ei = 62 ~ e and 63 = t — 2e it follows that 

gp{0,tM = \{l + lArcTH{t-2e)/t]^^ (4.12) 

for 2e < i < 3e, so we confirm the third equation in Eq.(4.2). 

It will be useful for the three projection case below to record the following related results. 
A similar analysis to that above yields 



r+_(ei,e2,e3) = / dyi [ dy2 — }. exp ( - 

Jo J-00 7r^eie2e3 V 



yl _ {yi - y2f _ yl_ 

Cl €2 ^3 



^ ' TT - 2ArcTan [J . ) ) (4.13) 



47r^(ei + €2 + €3) V VV ^2(61 + £2 + €3) 

In particular we then have 

r++ = r++(6,e,e) = -^ (4.14) 

3v3e 

We are able to evaluate the three projection case only in the situation where all time 
intervals are equal, so we set f — ne = e in the definition of Tc in Eq.(4.5). In this case, Tc 
possesses a number of helpful symmetries. The first is "refiection" symmetry: if we define 
—C by the string obtained by letting (H — >■—),( )■ +), (0 ^ 0) then 



-c — 



(4.17) 



The second symmetry is "time reversal" : if we define C by the string obtained by reversing 
the order of C, then we have 

Tc = Tc (4.18) 
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In addition we have the simple property that 



T.o... = + (4.19) 

These properties imply the following for the three projection case. On the one hand we 
have, 

= - = - + 7^+- (4.20) 

but we also have, 

r+++ = 7Vo+ - = r+o+ - (4.21) 
Combining these expressions gives, 

2r+++ = 7Vo+ + - To+_ (4.22) 

For each of the objects on the right, we may carry out the full range integral, to leave 

ourselves with an object of the form Tc{€i,e2,es) computed above. For example we have 
that T+OH- = ^++(e, 2e, e). Each of these objects may then be evaluated using Eq.(4.11) and 
Eq.(4.13), and combined to give, 

= -i= (4.23) 
4v4e 



gp{0,tM = ](^]' (4.24) 



It follows that 

nJn.t\nA))= , 

4 \27iit. 

when t = 4e, thus confirming the fourth equation in Eq.(4.2). 

To end this section, we confirm our claim in Section 3 that the function gp{0,t\0,0) drops 
to half its value immediately after a projection. On the face of it, this involves interpreting 
the expression (0|P which is ambiguous, so must be defined by a limiting procedure (where, 
recall (0| denotes {x\ at a; = 0). We thus consider the limit e„ — >■ in the expression 

gp{0, r|0, 0) = {0\e-'"'"Pe-'"'P ■ ■ ■ e-'"'Pe-'"'°\0) 

= J dy {0\e-'^'-Pe-'"'\y){y\P---e-'^'Pe-'"''\0) (4.25) 

We seek to show that the result obtained by this hmit is half the result obtained when the 
final projection is absent. We write 

P = ^{P + P) + ^{P-P) (4.26) 
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so that 

{Q\e-'"'^Pe-'^'\y) - \{Q\e-'^'-e-'^'\y) + \{Q\e-'^'"{P - P)e-'^'\y) (4.27) 

The first term on the right-hand side yields the sought after factor of a half as e„ 0, so it 
remains to show that the second term is zero. We have 

POO 

{0\e~'^'-{P - P)e-'"'\y) = dx (0|e-'^^"|a;) {{x\e-'^'\y) - {-x\e-'^'\y)) (4.28) 

Jo 

where we have used the fact that 

(0|e-^^^"|x) = (0|e-^^^"| - x) (4.29) 

(for the free particle, considered here). As e„ — >■ 0, (0|e~'^^"|a;) — >• 5{x) and the pair of terms 
in brackets on the right-hand side cancel in Eq.(4.28), so we obtain zero as required. 

V. ASYMPTOTIC LIMIT 

We now consider the asymptotic evaluation of 5'p(0,r|0,0) for large n in the case e„ = 
e = eg. This will give the values of the peaks of fp{t), Eq.(3.20). 

Recall that in the limit n — >■ oo and e ^ with t — {n+ l)e fixed, gp{xi, t\xo, 0) tends 
to the restricted propagator ^'^(^i) '''I^^O) 0). To obtain the asymptotic form of gp{0,T\0,0) 
for large n we therefore need to determine the lowest non-trivial correction to this result. It 
is clear from the definition Eq.(3.2) of gp that close to the hmit, we have the general form 

gp{xi,T\xo, 0) = gr{xi, t\xo, 0) + egi{xi, t\xo, 0) + O(e^) (5.1) 

for some function gi{xi,T\xQ, 0). Since the restricted propagator vanishes if either xi = or 
xq — 0, the object we need to calculate is 

j.(0,.|0,0)= Um ?£<^^ (6.2) 

A standard and convenient way to do this is to rotate to imaginary time f (we use tilde 
to denote Euclideanized time) and then write gp{xi,f\xo,0) as a Euclidean path integral. 
This is then defined in terms of the continuum limit of probabilities of random walks on 
a space time lattice of temporal spacing Af and spatial spacing rj. The details of this 
construction are very conveniently given by Hartle [23] so we will give only the briefest 
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account here. Using this language, g'p(0, f|0, 0) is then the continuum hmit of the object 
(277)~^iip(0, f|0, 0), where wp(0, f|0, 0) is the probabihty for a random walk on the lattice 
starting at the origin and ending at the origin at time f, with the restriction that the walker 
is in the positive x-axis at the intermediate times e, 2e, 3e ■ ■ ■ (where e > Af). 

The calculation of defined in this way, for values of e generally greater than the lattice 
spacing Af, is in fact a known problem in combinatorics called the tennis ball problem. It 
appears to have a formal solution, but this solution is too implicit for us to extract a useful 
result [24]. 

Fortunately, however, for the purposes of calculating the hmit Eq.(5.2) the results of 
Ref. [23] are sufficient. For this case, wc set e = Af and np(0, f |0, 0) is then the probability 
for a random walk from the origin to itself, with the restriction that the walker is in the 
positive X-axis at every intermediate step. On a finite lattice Hartle's calculations give the 
result, 

-MO,t|0,0)=(— ) . (5.3) 

to leading order for small e, rj [23]. We may use this result to compute the limit Eq.(5.2), 
which, continued back to real time, is 

9.0.-10.0) =(^) - (5.4) 
Through Eq.(5.1), this confirms Eq.(3.20) for large k. 

VI. A TIME AVERAGED RESULT 

We have argued that at the peak values of gp{0, t|0, 0) with n projections, we have the 
approximate result 

/ 77? \ 1/2 1 

(0|e-^^-P(n6)---P(26)P(6)|0)^^ (^) ^ (6.1) 

We showed above that this is exact for n = 1, 2, 3 and true asymptotically for large n. Some 
interesting exact results for any n may be obtained by considering a more general version 
of this object in which the projections are not restricted to act at the given set of times, 
e, 2e, 3e • • • . 

On the one hand, Eq.(3.17) may be expanded as a power series in powers of Vq: 

/ n1/2 ~ (—1)"' 
9.(0.r|0,0)= — E(bl^^oV (6.2) 

n=0 ^ ^ 
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On the other hand, the evolution operator with complex potential may be expanded in the 
usual S-matrix expansion, 

^-i{H-iV)r ^ e-'"^ T exp dtV{t)^ 

y^^-L / dtn... / dhT[V{tr,)...V{h)] (6.3) 

rr^ ^- Jo Jo 



-iHt 



e 

n=0 



where T denotes time ordering, V{t) = VoP{t) and P{t) = e'^^Pe'^^K It follows that 

^y(0,r|0,0) = (0|exp(-^(i/-^\/)r)|0) 

= y^^^V^ / dtn... / dh (0|e-^^^r[P(y...P(ii)]|0) (6.4) 

n=0 

Equating powers of Vq in Eqs.(6.2) and (6.4) and writing out the time ordering explicitly, 
we deduce that 

I'tn r^2 



- / dtn dtn-1... / dh (0|e-^^^P(i„)...P(ii)|0) = (-— ) ^ (6.5) 



We would like to write this in a form involving P, instead of P. We introduce the reflection 
operator 

(6.6) 

and note that P = RPR, the Hamiltonian H commutes with R, and P|0) = |0) (where, 
recall, |0) denotes \x) at x = 0). It follows that 

(0|e-^^-P(tJ...P(ti)|0) = (0|e-'^-P(O...P(ti)|0) (6.7) 

so that Eq.(6.5) holds with all the P's replaced with P's. Noting that 

— / dtn dtn-1... / dh^l (6.8) 

T Jo Jo Jo 

we see that Eq.(6.5) is of the desired general form, Eq.(6.1), but time-averaged over the 
times of the projections. 

The question is now to what extent the time-averaged expression on the left-hand side of 
Eq.(6.5) is close to Eq.(6.1). We may take this further in two different ways. 

First, note that for a real-valued function / of n variables, we have the mean value 
theorem 

— / dtn dtn-i- / dtif{tn, t^-i ■ ■ ■ h) ^ /(^„, • • • ^i) (6.9) 
Jo Jo Jo 
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for some set of numbers ^„ > > • • • > in the interval [0, r]. The integral in Eq.(6.5) 
is easily made into a real integral over a real- valued function by analytic continuation, and 
we therefore deduce the exact result 

(0|e— P«„)...P(6)|0>=(^)"''^^ (6.10) 

for some set of projection times ^„ > > ■ ■ ■ > ii 'va. the interval [0, r]. 

Second, we may expand the integrand in the left-hand side of Eq.(6.5) about the values 
ke, to get some insight into why these particular values have any special significance. We 
have 

(0|e-'^^P(i„)...P(ii)|0) = (0|e-^^^P(ne)...P(e)|0) 

k=l ^ 

+••• (6.11) 
Inserting in Eq.(6.5), and noting that 

dtn / dtn-1--- / dti tk = ke (6-12) 
JO Jo 

(where, recall, r = (n-|-l)e) we see the first order term in the expansion Eq.(6.11) averages to 
zero in Eq.(6.5). This is clearly only true for expansion about the special values tk = ke. We 
therefore obtain the desired result Eq.(6.1) up to second order corrections. This suggests 
that the values tk = ke are significant because they give the best approximation to the 
average in Eq.(6.5). These results give evidence that Eq.(3.20) holds approximately for all 
n, including the intermediate values not covered in the previous two sections. 




VII. NUMERICAL RESULTS 

To support the analytic results described in the previous sections we evaluate Eq.(3.18) 
numerically for up to 20 projections and confirm the conjectured form Eq.(3.19), (3.21). It 
is convenient to define a sequence of functions Fn{t,x) defined by 

Foit,x) = {x\e-''''\0) (7.1) 

for < t < e and 

Fnit, x) = (x|e-'^(*-"^) (Pe-^^^)" |0) (7.2) 
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for ne < t < {n+ l)e where n — 1, 2, 3 • • • . The desired object Eq.(3.18) (with, for conve- 
nience, eo = e), is then given by 

5p(0,i|0,0) = F„(i,0) (7.3) 
The sequence Fn{t,x) may be calculated using the recursion relation, 

/■oo 

F^{t,x)^ / (x|e-^(*— )|y)F„_i(ne,?/) (7.4) 
Jo 

Using a new time coordinate s defined hy t — se, rotating to imaginary time, and defining 
Fn{s,x) — Fn{—it,x), we have 

Fn{s,x) = ( ' Hdy exp ( -!!$^^] F^_^{n,y) (7.5) 



2n{s — n)J Jo V 2(s — n) 

for n < s < n + 1. 

We have attempted to find an approximate analytic solution to Eq.(7.5) for large n, 
but without success. However, a numerical solution is straightforward and yields all the 
information we require. This was done using a simple mid-point rule. The lattice size was 
chosen to be of order 10~^, and the results were checked for robustness against changes in 

lattice size. 

The numerical result for fp{t) (defined in Eq.(3.19) in terms of gp{0,t\0,0)) is plotted 
in Figure 3, along with our claimed approximate analytic expression for fp{t), Eq.(3.21). 
We see that there is excellent agreement. The values at the peaks and troughs appear to 
agree perfectly. The only small discrepancy is that the interpolating functions between the 
peaks and troughs are not exactly linear. This discrepancy is only noticeable for intermediate 
values of n and in any event since, as argued, the curve is effectively averaged over time in the 
PDX, this discrepancy is insignificant. We therefore find substantial numerical confirmation 
for our our main result, Eqs.(3.19), (3.21). 

The apparently perfect agreement of numerical results with the approximate analytic 
expression Eq.(3.21) at the peaks and troughs is striking. We wonder if the approximate 
analytic expression is in fact exact at these points, but we have not been able to prove this, 
except for the cases n = 1, 2, 3 and for large n. 

A useful way of seeing even more precisely the relationship between fp{t) and fv{t) is to 
define the function 

Jv[t) 
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so that 

fp(t)^{l + S(t))Mt) (7.7) 

This is plotted in Figure 4. It is a simple function oscillating around zero between ±1/3 
with period e. In terms of S{t), the relationship between gp and gv then has the particularly 
simple form 

gp(0, t\0, 0) = (1 + S(t))gv(0, t\0, 0) (7.8) 

This relationship is perhaps the most concise summary of the sought-after connection be- 
tween the propagators. 



VIII. TIMESCALES 



We now give a more detailed explanation as to the timescales involved in proving the 
approximate equivalence of gy and gp. We broadly expect that the appropriate timescale 
is the Zeno time of the initial state, tz = 1/AH. However, we have derived a very precise 
connection between evolution in the presence of a complex potential and evolution with 
projection operators so we are in a position to investigate the specific way in which an 
initial state may discriminate between these two types of evolution. 

We have argued that the equivalence boils down to proving the equivalence of propaga- 
tion along the boundary, Eq.(3.12), and we have shown that gp{0,t\0,0) oscillates around 
gv{0,t\0,0) with period e. Suppose we have an initial state iIj{x,0). Let us consider the 
change in the wave function arising from replacing gv with gp using the PDX, Eq.(3.10). It 
is 

6ilj{x,,r) = dt2j^'dti ^{x^,T\xM)\,=^5g{QM^M) |^(0'^i) (8-1) 

where 

5^(0, i2|0, h) - (0, i2|0, h) - ^p(0, t2% h) (8.2) 

and we also used Eq.(3.11). The important part of this expression is the ti integral, which 
we expect will be small if the initial state is sufficiently slowly varying in time. 

The results of Section 3 and 8 show that dg{0,t\0,0) oscillates around zero with period 
e. In particular, Eq.(7.8) shows that 

6g{0,t\0,0) = S{t)gv{0,t\0,0) (8.3) 
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The explicit form of S{t) is given in Figure 4, but its important qualitative feature is its 
oscillation with period e, so for simplicity, S{t) may be loosely modelled by the function 
sin(27rt/e). To be definite we take the initial state to be a Gaussian wave packet, so we 

have 

^|;{x,t) = Nexp(^-^^^^^^^^^^^ + ipx-iE?j (8.4) 

where E — p^/2m, N is a. normalization factor, and we have ignored wave packet spreading 
effects in evolving the state. For such a state the Zeno time is tz = ma /p. In Eq.(8.1) 

the differentiation of ip produces a prefactor which does not contribute to the leading order 
evaluation of the time integral, and if the range of integration is much greater than e, we 
obtain the order of magnitude result 

dh 6g{0,t2\0,t,) 5^(0, ti) - exp ( -%{Ee - lA (8.5) 



If -Ee <^ 1, the right-hand side is clearly very small if e <^ t^. If i?e > 2, the right-hand side 
if bounded from above by exp(— t|/e^) so again will be small if e ^ tz- Hence the Zeno time 
of the initial state is indeed the timescale controlling the validity of the approximation, as 
expected. Note also that Eq.(8.5) goes to zero as e — > 0, as it must, since gp and gv become 
exactly equal (and equal to gr) in this limit. 

The only problematic case is that in which the initial state has energy £■ ~ 1/e. In this 
case, the right-hand side of Eq.(8.5) is not necessarily small and the approximation may fail. 
This is not surprising since it is the case in which the oscillations in time of the incoming 
state are comparable to the temporal spacing of the projections, so that the state can "see" 
the difference between the complex potential and projections at a discrete set of times. 

Hence, apart from the above exception, Eq.(3.12) holds for e <^ tz- As outlined in the 
Introduction, reflection in the complex potential is negligible if Vq <^ £^ which is equivalent 
to e ^ l/E". Therefore there is an interesting regime, namely 

I « e « ^ (8.6) 

in which the approximate equivalence Eq.(3.12) holds, yet there is neghgible reflection. 
This regime is important in, for example, study of the arrival time problem using complex 
potentials [9]. This issue will considered in more detail elsewhere [26]. (See also Ref.[25] for 
an interesting discussion of timescales in the Zeno effect.) 
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IX. SUMMARY AND DISCUSSION 



This paper was physically motivated by the desire to understand the effect of periodically 
acting projections onto the positive x-axis for a free particle, Eq.(l.l). A valuable way to 
proceed is to use the conjectured relationship Eq.(1.3) with a complex potential first put 
forward by Echanobe et al. This connection, together with known results on scattering, 
establishes the timcscalc under which significant reflection occurs in Eq.(l.l). We noted 
that the arguments for the relationship Eq.(1.3) are only heuristic and there is scope for a 
more substantial proof. 

We proved Eq.(1.3) by considering the associated propagators Eqs.(3.1), (3.2). We noted 
that an approximate equivalence between these propagators is expected since both prop- 
agators tend to the restricted propagator gr{xi,T\xo,0) in the limits e ^ 0, n ^ oo and 
Vq ^ oo. The path decomposition expansion reduced the proof of equivalence of these prop- 
agators to the simpler case of propagation between points lying on the origin, Eq.(3.12). 
The propagator along the origin for the complex potential gv{0,t\0,0) is already known, 
Eq.(3.17), so the bulk of the proof was to derive the analogous result for the propagator 
with projections, gp{0,T\0,0). Our main result was to prove that this propagator has the 
approximate form Eqs.(3.19), (3.21), which we proved using a variety of analytic and nu- 
merical methods. In effect, the main achievement of this paper has therefore been to obtain 
a good approximate analytic expression for the propagator gp that appears in Eq.(l.l). 

We found that gp{0,t\0,0) oscillates with period e around gv{0,t\0,0) as long as Voe 
4/3, a result most concisely summarized in Eq.(7.8). The approximate equivalence Eq.(3.12) 
of these propagators then holds in a time-averaged sense as long as the timescale e between 
projections is much smaller than the Zeno time of the initial state, 1/AH (but may fail in 
the special case when the incoming state is peaked about energy E ~ 1/e)- These conditions 
agree in essence with those of an extended version of the results of Echanobe ct al, with an 
advantage over their results that a definite relationship between Vq and e is obtained. We 
noted that their restriction Eq.(1.4) relating Vq and e is in fact stronger than required and 
the equality Voe 4/3 derived here gives the best approximate equivalence between gp and 

gv- 

In addition to helping establish the timescales for reflection in Eq.(l.l), the connection 
Eq.(1.3) has been of use in investigating the near-Zeno regime of Eq.(l.l) [14-16]. In the 
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limit e — > 0, the Zeno limit, motion in Eq.(l.l) is confined entirely to states with support 
only in X > and is described by the restricted propagator Eq.(3.8). It is of interest to 
explore the nature of the dynamics very close to this regime, i.e., for small but finite e. 
This can be achieved by examining the form of the propagator with complex potential gy 
for large Vq. In this regime, one can calculate, for example, the (unnormalized) probability 
density for crossing the origin between r and t + dr, which is found to be 

dN 



n(r) 



dr 

^ {i;fiT)\p5{x)p\Mr)) (9.1) 



Here, N{t) is the survival probability (given by the norm of the state evolved in the presence 
of the complex potential) and |'0/(t)) is the freely evolved wave function. This is proportional 
to the average kinetic energy density at the origin. As anticipated it goes to zero as Vq — > oo 
since there is total refiection. However, it is useful to define a normalized distribution, 

n (.^ - "(^) 

Jo dtU{t) 

where (p) is the average momentum in the initial wave packet [14-16]. Interestingly, this is 
now independent of the complex potential despite this being the regime of strongly-acting 
measurement. This regime will be explore further in future publications [26]. 

Finally, we comment on the possible generality of the connection Eq.(1.3). We first note 
that it may in fact be written 

where P — 1 — P. We have proved Eq.(9.3) for the case in which the projections are onto the 
positive X-axis, but it seems reasonably clear that the relationship will hold for projections 
onto any region [a, b] of the x-axis (as long as it is not too small) and indeed for regions of 
configuration space in a many-dimcnsional model. Such a potential has been used recently in 
the decoherent histories analysis of quantum cosmological models [27]). Moreover, although 
the proof of Eq.(9.3) given in this paper relied heavily on the fact that P projects onto 
position, the form Eq.(9.3) and the heuristic argument for it in Section 2 do not rely on the 
particular form of P. We therefore conjecture that Eq.(9.3) may hold for a wider variety of 
projection operators, not just projectors onto position. This will be pursued elsewhere. 
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Figure 1. The path decomposition expansion for gp. Any path from > at i = to a 
final point > at t = r which crosses a; = has a first crossing of x = at ti and a last 
crossing at ^2- The crossing part of the propagator from (a;o,0) to (a;i,r) may therefore be 
decomposed into three parts: (A) restricted propagation entirely in a; > 0, (B) propagation 
starting and ending on x = with the restriction that the paths are in x > at a discrete 
set of times, and (C) restricted propagation entirely in x > 0. The corresponding path 
decomposition expansion formula is given in Eq.(3.10). (The Qr term corresponds to paths 
which never cross a; = 0). 
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Figure 2. A plot of the functions fv{t) and fp{t), defined for the complex potential by 
Eq.(3.17) with V^e = 4/3, and for intermittent projectors by Eqs.(3.19), (3.21). (The time 
scaling is such that e = 1 in the plot, so that the peaks occur at integer values of t). We 
see that fp{t) oscillates around fv{t) with period e so the two functions are equal in a time- 
averaged sense, when integrated against functions which are slowly varying on a timescale 
e. 
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Figure 3. A plot of the numerical calculation of fp{t) (the set of dotted lines), and the 
conjectured analytic expression for fp{t), Eqs.(3.19), (3.21) (the set of straight lines). (Again 
we use a time scaling such that e = 1). There is excellent agreement between analytic and 
numerical results at the peaks and troughs of the functions. The numerical result shows 
some departure from exact linear behaviour between the peaks and troughs but this is 
insignificant as argued in the main text. A plot of fv{t), the smooth curve, is also given for 
reference. 
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Figure 4. A plot of the function S{t) defined in Eq.(7.6), describing the oscillations of 
5'p(0,t|0,0) around gv{0,t\0,0). The bold line represents the numerical calculation and the 
dashed line the analytic result (from Eqs.(3.17), (3.21)). After the first few oscillations, it 
oscillates between ±1/3 (approximately). 
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